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Abstract — We consider a networked control system where 
a linear time-invariant (LTI) plant, subject to a stochastic 
disturbance, is controlled over a communication channel with 
colored noise and a signal-to-noise ratio (SNR) constraint. 
The controller is based on output feedback and consists of 
an encoder that measures the plant output and transmits over 
the channel, and a decoder that receives the channel output 
and issues the control signal. The objective is to stabilize the 
plant and minimize a quadratic cost function, subject to the 
SNR constraint. 

It is shown that optimal LTI controllers can be obtained by 
solving a convex optimization problem in the Youla parameter 
and performing a spectral factorization. The functional to 
minimize is a sum of two terms: the first is the cost in the 
classical linear quadratic control problem and the second is a 
new term that is induced by the channel noise. 

A necessary and sufficient condition on the SNR for stabi- 
lization by an LTI controller follows directly from a constraint 
of the optimization problem. It is shown how the minimization 
can be approximated by a semidefinite program. The solution 
is finally illustrated by a numerical example. 

Index Terms — ACGN channel, control over noisy channels, 
linear-quadratic-Gaussian control, networked control systems, 
signal-to-noise ratio (SNR) 



I. Introduction 

COMMUNICATION limitations are a fundamental 
characteristic of networked control. The recent trend 
of decentralized and large-scale systems has therefore 
driven an interest in research on how control systems are 
affected by communication phenomena such as random 
time delays, packet losses, quantization and noise. A 
popular approach used for research on the fundamental 
aspects of communication limitations in control systems is 
to consider a plant that is controlled over a communication 
channel, as depicted in Fig. [TJ Due to the lack of a theory 
that can handle all limitations at once, the channel model 
is typically chosen to highlight a specific issue. 
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Fig. 1. Feedback control of a plant, with disturbance, over a 
communication channel. The control system consists of an encoder, 
which also does measurement filtering, and a decoder, which also 
determines the control signal. 



One example is given by data-rate constraints, whose 
study has led to one of the most well-known results in the 
area, known as the data-rate theorem. In discrete time, 
it says that an unstable linear plant G can be stabilized 
over a digital error-free channel if and only if 



K > Vmax{0,log 2 |A,(G)|} = H G , 






(1) 



where 7Z is the rate of the channel and Xi(G) is the ith 
pole of G [J, 0, 0. 

The situation is more complicated for noisy channels. A 
necessary and sufficient condition for almost sure asymp- 
totic stabilizability is that the channel capacity C satisfies 
C > 'Hg [!]• But since this condition is not generally 
sufficient for mean-square stability, the concept of any- 
time capacity has been proposed for characterization of 
moment stabilizability [5J. 

Stability is, however, easier to characterize for control of 
a linear time-invariant (LTI) plant over an additive white 
Gaussian noise (AWGN) channel or, more generally, an 
additive colored Gaussian noise (ACGN) channel. Since 
the communication aspect highlighted by this channel 
model is a Signal-to-Noise Ratio (SNR) constraint, this 
setting will be referred to as the SNR framework. 

The SNR framework is mainly attractive due to its 
simplicity. However, it has been argued that the usage 
of linear controllers admits application of established 
performance and robustness tools [6_ . The obtained results 
can sometimes also be used to draw conclusions about 
and design controllers for other communication limitations 
such as packet drops [7] or rate limitations [5], [9]. More- 
over, the SNR framework can be useful for applications 
such as power control in mobile communication systems 
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A. Main Result and Outline 

This paper considers the problem of control design 
in the SNR framework. The system has the structure 
illustrated in Fig. Q] The plant is LTI, possibly unstable 
and subject to a Gaussian disturbance. The controller is 
based on output feedback and consists of an encoder and 
a decoder. The encoder measures the plant output, filters 
the measurements and encodes them for transmission 
over an ACGN channel. The decoder receives the channel 
output, decodes it and determines the control signal. The 
objective of the controller is to stabilize the system and 
minimize a quadratic cost function, while satisfying the 
SNR constraint. 

The main result is that an optimal LTI output feed- 
back controller can be obtained by minimizing a convex 
functional and performing a spectral factorization. The 
minimization is performed over the Youla parametrization 
of the product of the encoder and the decoder, and the 
functional is the sum of the classical LQG cost and a new 
term that is induced by the channel noise. A condition for 
stabilizability, which coincides with the previously known 
condition in the AWGN case, is obtained as a constraint 
of the minimization problem. It is shown how to formulate 
the minimization as a semidefinite program. As a by- 
product of the main result, it is also shown how the 
encoder and decoder should be chosen in order to minimize 
the impact of the channel noise while preserving the closed 
loop transfer function given by a nominal LTI controller 
that has been designed for a classical feedback system. 

The rest of this section will present the previous research 
in the SNR framework. Section [TT] presents the mathe- 
matical notation used in this paper. The exact problem 
formulation is given in Section Mil Section IIVI is devoted 
to the solution of the problem. Section |V] presents a 
procedure for numerical solution and a numerical example. 
Finally, Section |VT1 presents the conclusions and discusses 
further research. Some technical lemmas have been put in 
the appendix. 

B. Previous Research 

Necessary and sufficient conditions for stabilizability, 
similar to the data-rate theorem, have been found for the 
SNR framework. They do, however, vary depending on 
some of the assumptions. Generally, the condition for the 
AWGN channel is that the SNR a 1 satisfies the inequality 



a 2 > 



n 



|max{l,Ai(G)}| 2 ) -l + V + S, 



(2) 



where r\ and S depend on the specific assumptions, as will 
be explained. Assuming no plant disturbance and static 
state feedback, the condition is that (J2| holds with r\ = S = 
0. By writing fl} and @ in terms of the respective channel 
capacities, it can be shown that the capacity requirements 
for stabilization in the two settings are equal [BJ. 



For LTI output feedback, the condition is again @ but 
now the terms r\ and <5 are non-negative and depend on 
the non-minimum phase zeros and the relative degree of 
the plant, respectively [BJ. If there is a plant disturbance, 
the same condition holds if the controller has two degrees 
of freedom (DOF) [11] but not if it only has one DOF, in 
which case the required SNR is larger [12]. The condition 
with T] = S = can be recovered for the output feedback 
case, either by introducing channel feedback, meaning that 
the encoder knows the channel output |llj , or by allowing 
a time- varying controller, although the latter leads to poor 
robustness and sensitivity |6J. Further, it has been shown 
that the condition ([2]) with r\ = S = is necessary for 
stabilizability even if nonlinear and time varying state 
feedback controllers are allowed |13| . 

Similar conditions have been found for LTI control of 
a plant with no disturbance over an ACGN channel [12J. 
The case with first order moving average channel noise 
was further analyzed in [14] . 

An early formulation of a feedback control problem 
over an AWGN channel with feedback was made in [15] . 
It was shown how to find the optimal linear controller 
and that it is globally optimal for first-order plants. A 
counterexample was provided, showing that non-linear 
solutions may outperform linear ones for higher order 
plants [T3]. It should, however, be noted that the provided 
counterexample requires the plant to be time- varying and 
that the encoder has a memory structure where it does 
not remember past plant output. 

Since then, many authors have considered similar con- 
trol design problems that have been simplified by as- 
sumption of a certain controller structure, see |16| . [17] , 
[18j . |19) . |20| . A quite general approach was proposed 
in [11] . but it was also noted that it leads to a difficult 
optimization problem with sparsity constraints when it is 
applied to controllers with two degrees of freedom. 

The problem of optimizing the control performance at 
a given terminal time was considered in [21] and [22]. The 
solutions may however yield poor transient performance 
and can therefore be unsuitable for closed-loop control. 

A lower bound on the variance of the plant state was 
obtained for feedback control over AWGN channels, using 
general controllers with two degrees of freedom, in |13| . 
This bound tends to infinity as the SNR approaches the 
limit for when stabilization is possible. 

An important contribution was recently made in [9]. 
Although the paper mainly considers control over a 
rate-limited channel, this is done through design of an 
LTI output feedback controller in the SNR framework. 
assuming an AWGN channel with feedback. The optimal 
performance is shown to be obtained by solving a convex 
optimization problem with the same structure as the 
one obtained in this paper. An optimal controller is 
then acquired by finding rational transfer functions that 
approximate certain frequency responses. Related results, 
for when the controller is pre-designed and the coding 
system should have unity transfer function, are given in 
[5] and [23]. 



Comparing with the results presented here, the case 
without channel feedback is not mentioned in [SJ. The 
presented convex functional that gives the optimal cost 
for the case with channel feedback can, however, be 
modified to give the optimal cost for this case as well. 
The expressions for the optimal transfer functions that 
are given can, with additional work, also be modified to 
give solutions to the case without channel feedback. 

We claim that the solution presented in this paper has a 
clearer structure than the one given in [9] . For example, we 
do not require an over-parametrization of the controller. 
Moreover, while the plant is assumed to be single- input 
single-output (SISO) in [SJ, it is here allowed to be slightly 
more general, making it possible to include any number 
of noise and reference signals and to penalize the control 
signal variance. Also, we allow the channel noise to be 
colored. A final contribution of this paper relative to [pj 
is that it is shown how to pose the optimization problem 
as a semidefinite program. 

The approach used in this paper is based on the solution 
of a communication-theoretic problem involving design 
of encoders and decoders for a Gaussian source and a 
Gaussian channel when there is a delay constraint |24j . 
Some instances of that problem can be viewed as special 
cases (open loop versions) of the problem considered here, 
which therefore may be viewed as a partial generalization. 

II. Notation 

The proofs in this paper make extensive use of concepts 
from functional analysis, such as C p (Lebesgue), T-L p 
(Hardy) and N + (Smirnov) function classes and inner- 
outer factorizations. To conserve space, only some of the 
most important facts will be given here. The interested 
reader is referred either to [25] or to [26], [27] and [28] for 
the remaining relevant definitions and theorems. 

The natural logarithm is denoted log. The complex 
unit circle is denoted by T. For matrix-valued functions 
X(z),Y(z) defined on T, define 
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where \\-\\ F is the Frobenius norm and o\ the largest 
singular value. 

A transfer matrix X(z) is said to be proper if the 
mapping z t— > X(l/z) is analytic at 0. It is strictly proper 
if also lim z ^ooA(z) = 0. The space of all rational and 
proper transfer matrices with real coefficients is denoted 

by n. 

Cp, for p = 1,2,00, is the space of matrix- valued 
functions X(z), defined on T, that satisfy ||A|| < oo. 



G 



C 



-K+ 



D 



H 



ii- 



Fig. 2. Feedback system with ACGN communication channel. The 
objective is to design the controller components C and D so that 
the plant G is stabilized and the variance of z is minimized under 
the SNR constraint E(£ 2 ) < a 2 . H is a spectral factor of the channel 



The subspaces 1Z£ P consists of all real, rational and proper 
transfer matrices with no poles on T. 

H pi for p = l,2,oo, is the space of matrix-valued 
functions X(z) such that z h-> X(1/z) is analytic on the 
open unit disk and 

s up||A r || < oo, 

■r>l 

where X r (z) — X(rz). The subspaces 7ZH P consists of all 
real, rational, stable and proper transfer matrices. When 
a function in Up is evaluated on T, it is to be understood 
as the radial limit lim r _ > . 1 + X(rz). 

The arguments of transfer matrices will often be omitted 
when they are clear from the context. Equalities and 
inequalities involving functions evaluated on T are to be 
interpreted as holding almost everywhere on T. 

III. Problem Formulation 

A detailed block diagram representation of the system is 
shown in Fig. [2] The plant G is a multi-input multi-output 
(MIMO) LTI system with state space realization 



G(z) 



G zv (z) G zu (z) 

Gy V [Zj Gy U [Z) 
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C 2 



B\ B 2 



Du D 12 
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where (A, B 2 ) is stabilizable and (C 2 ,A) is detectable. 
The signals v and z are vector-valued with n v and n z 
elements, respectively. All other signals are scalar- valued. 
Accordingly, G zv is n z x n v , G yv is 1 x n v , G zu is n z x 1 
and G yu is scalar and strictly proper. It is assumed that 
G* ZU G ZU and G yv G* v have no zeros or poles on T. 

The input v is used to model exogenous signals such 
as load disturbances, measurement noise and reference 
signals. It is assumed that v and w are mutually indepen- 
dent white noise sequences with zero mean and identity 
variance. The other signals in Fig. [2] are the channel noise 
n, the control signal u, the measurement y and the control 
error or performance index z. 

The feedback system is said to be internally stable if 
no additive injection of a finite-variance stochastic signal 



*0 



G ,, 



II 



c 



D 



at any point in the block diagram leads to another signal w\ 

having unbounded variance. This is true if and only if all 
closed loop transfer functions are in 'H.i . 

The communication channel is an ACGINLj channel with 
SNR er 2 > 0. The channel noise has the spectral factor 
H £ TZHoo, which is assumed to have no zeros on T. 
Since the channel input and output can be scaled by C L 

and D, it can be assumed that n has unit variance and thus 

\\H\\l = 1 without loss of generality. The SNR constraint Fig 3 Block diagram for internal stab iiity analysis, 
is then equivalent to the power constraint 

E(t 2 ) < a 2 . 



W-2 



(3) 



The objective is to find causal LTI systems C and D that 
make the system internally stable, satisfy the constraint 
(|3|) and minimize E(z T z) in stationarity. 

By expressing z and t in terms of the transfer functions 
in Fig. [2l the objective and the SNR constraint can be 
written as 



J(C,D) 



and 



G, 
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DCG ZU G V 



1 - DCG V 



DUG, 



1 - DCG 
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DCHG 
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1 - DCG 



1 - DCG,, 



<a A 



yu. 



(4) 



respectively. The main problem of this paper is thus to 
minimize J(C, D) over C and D subject to (j4]) and internal 
stability of the feedback system. 

For technical reasons, only solutions where the product 
DC is a rational transfer function will be considered. This 
may exclude the possibility of achieving the minimum 
value, but the infimum can still be arbitrarily well approxi- 
mated by such functions. Since D and C are required to be 
proper, DC has to be proper as well. That is, DC £ 71. 
Though the latter will be enforced, it is not explicitly 
required that C and D are proper. It will, however, be 
seen that the solution is constructed so that C £ T-L 2 is 
outer. Then C, C _1 are proper, and D = (DC)C~ 1 is also 
proper. 

IV. Optimal Linear Control 

The solution of the problem presented in the previ- 
ous section is divided into three subsections. The first 
characterizes internal stability of the system. The second 
introduces the optimal factorization of a given nominal 
controller. The third section shows that the optimal 
factorization result can be used to find an equivalence 
between the main problem and the minimization of a 
convex functional in the Youla parameter. 

A. Internal Stability 

The product DC will play an important role in the 
solution. Therefore, introduce 

K = DC. 

1 Since only linear controllers are considered, it does not matter if 
n or v are Gaussian or not. Linear solutions may, of course, be more 
or less suboptimal depending on their distributions. 



Following the same reasoning as in [29], it is concluded 
that internal stability of the systems in Fig. [2] and Fig. [3] 
are equivalent (H does not have to be included since it 
is open-loop stable and not part of the feedback loop). 
The latter can be represented by the closed loop map T, 
defined by 



v 

t 


= T 


Itfl 
U>2 


a. 




n 



Hence, the system in Fig. [5] is internally stable if and only 
if 



T = 



KGy U 



DGy U 



l-KG v 
C 

I -KG 
K 



1 - KG VU 1 - KG VU 1 - KG % 



U" 



1 - KG yu 

CGy U 
1 - KGy U I 

KG VU 



1 — KGy U 
KGy U 

KGyu 

D 



en 2 



(5) 



The following two lemmas will give necessary and 
sufficient conditions for internal stability, respectively. 



Lemma 1: Suppose that T £ T-L 2 , G yu 



NM~ 



coprimc factorization over TZT-Loo and that U, V £ TZTi 
satisfy the Bezout identity VM + UN = 1. Then 

MQ-U 



K 



Q£KH C 



NQ + V 
Proof: It follows directly from ([5]) that 



G,, 



1-KGu 



£U 2 
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1 ~ KGyu 
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1 ~ KGyu 



(G) 



£U 2 - 



'yu x - li, -'yu x Ji, ~'yL 

These transfer functions are rational and have no poles 
on or outside the unit circle, so it follows that 



f 



-K 
1 



1 



K 



1 - KGyu 1 - KG. 



G 



yu 



1 



11 » 



1 — KGyu 1 — KGyu 



£KH C 



(7) 

The set of K satisfying ([7]) can be parametrized using the 
Youla parametrization of all stabilizing controllers [29|. 
That is, K can be written as in §Q. ■ 

Lemma 2: Suppose that 

MQ-U 



K = DC 



NQ + V 



Q£KH C 



(8) 



where G yu = NM 1 is a coprime factorization over TZU-c 
and U,V £ WH^ satisfy the Bezout identity VM + UN ■- 
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Fig. 4. Classical feedback system without communication channel. 



1. Suppose also that C € Hi is outer and that D £ £ 2 - 
Then T £H 2 . 

Proof: It follows from (El that 



G„ 



1 - KG yu 

KLryu 



£ KUr 



K 



1 - KGy U 



eWHc 



1 



1 — KGy U 1 _ KGyu 



leKHr- 



Moreover, 



DG 



yy. 



KG 



yu 



1 - KGyu 1 - KG 



-c- 



yv 



where the left hand side is in £2 since it is the product 
of an £2 function and a IZJioa function. Since C is outer, 
application of Lemma [7] (in the appendix) gives that the 
right hand side is in % 2. A similar argument shows that 



D 



I — KGyu 



£% 2 



Finally, 



C 



en 2 



LjGy U 
1 - KGyu 



€H 2 



1 — KGyu 

since these functions are products of an 14.2 function and 
an TZHoo function. Since TZHqo C "H 2 it has been proved 
that all elements of T are in H2 and so T £ %2- ■ 

B. Optimal Factorization 

Suppose for now that K £ 1Z is a given stabilizing 
controller for the classical feedback system in Fig.0] Thus, 
K satisfies ([6]). Nothing else is assumed about the design 
of K. It could for example be the I-L2 optimal controller 
or have some other desirable properties in terms of step 
responses or closed loop sensitivity 

In either case, it is a natural question to ask what the 
best way is to implement this controller in the architecture 
of Fig. [2] If the nominal design is to be preserved then 
C and D should satisfy K — DC since the transfer 
matrix from v to z would then be the same. Given 
this relationship, choosing C and D can be thought of 
as factorizing K. The factorization should be chosen 
to minimize the negative effect of the communication 
channel. That is, they should keep the system stable, 
satisfy the SNR constraint and minimize the impact of 
the channel noise. That is, to minimize the contribution 
of n to E(z T z). 



Rewriting J(C, D) and the SNR constraint g]) with DC 
replaced by K gives 



J(C, D) 



and 



G zv 



CGyy 



■** £* zu£*yv 



1 — KGyu 



DHG Z 



1 — K Gjm 



KHG yu 



I -KG 



1 — KGy U 



<a' 



11 » 



(9) 



(10) 



The objective of the optimal factorization problem is to 
find C and D such that ^ is minimized subject to (JTUJ) 
and K = DC. Stability is not considered now but it will 
be seen that the obtained solution is stabilizing anyway. 
The optimal factorization will later be used to solve the 
main problem of this paper. Alternatively, it could also be 
used to factorize a nominal K that was designed for the 
classical feedback architecture. 

Note that, for given K, the first term in © is constant 
and that the second term is a weighted norm of D. In 
the left hand side of (ITU)) , the first term is a weighted 
norm of C and the second is constant. Thus, the optimal 
factorization problem is a minimization of a weighted norm 
of .D, subject to an upper bound on a weighted norm of 
C and the constraint K = DC. 

Before the solution to this problem is given, it is noted 
that the SNR constraint will be impossible to satisfy unless 
K satisfies 



dcf O 

a = a 



KHG n 



l-KGu 



>0. 



Actually, if a = then, since G yv G* has no poles or zeros 
on T, 



CG™ 



1 — DCGy U 







C 







KHGyu 



0, 



1 — DGGyu 1 — KGyu 

which is a contradiction. Thus, it will be assumed that K 
is such that a > 0. Introducing 



1 



1-KGy 



£KH a 



for notational convenience, the set of feasible pairs (C, D), 
parametrized by K , is defined as 



9 



C,D 



(K) = {(C, 



D) 



\CSG 



yv\\ 2 



<a,DC = A'j 



The solution to the optimal factorization problem is now 
given by the following lemma. 

Lemma 3 (Optimal Factorization): Suppose a > 0, S £ 
KHoc, K £ 11 and that H £ TZHao, G* ZU G ZU £ TIC^ and 
G yv G* £ TZCoo have no zeros on T. Then 



inf \\DSHG 



1 2 

ZM||2 — II J^G ZU Gy V L 



(11) 

Suppose furthermore that K £ 1ZC\ satisfies (j6]). Then 
there exists (C,D) £ &c,d(K) with C £ H2 outer and 
D £ £ 2 , such that the minimum is attained and ([TT]) holds 
with equality. 



If K is not identically zero, then (G, D) is optimal if 
and only if DC = K and 



icr 






|"" "trzutrj/ulli y Gy V G yv 



KH\ on T. (12) 



If K = 0, then the minimum is achieved by D — and 
any G that satisfies ||C5G a „|| 2 < a. 

Proof: If if = then the right hand side of QTQ is 0. 
Letting D = gives ||DS'iIG zu || 2 = and it is clear that 
(C, -D) € ©c,-D if C is as stated. 

Thus, it can now be assumed that K is not identically 
zero. Then G is not identically zero and D = KC^ 1 . 

By assumption both G* ZU G ZU and G yv G* are positive 
on the unit circle. Since these functions are rational this 
implies that 



3e > such that G* ZU G ZU > e and G yv G 



yv 



> e, on T. 



(13) 

Thus by the factorization theorem in [30J there exist scalar 
minimum phase transfer functions G zu , G yv £ H2 such 
that 

^zu^zu — ( J r 2M ( - r z« 1 ^yv^yv — l ~ r yv'^ yv - 

Now, yC^Gy^ || 2 < a and Cauchy-Schwarz's inequality 
gives 

2 



\DSHG ZU \ 



> 



KG- X SHG ZU 
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KC^SHGzu 



Kb HG zu Gy V 
Kb HG zu G yv 



This proves the lower bound (|lll) . 

Equality holds if and only if \KC- l bHG zu \ 
\CSG yv \ are proportional on the unit circle 



yv\\ 2 



and 

and 

a. It is easily verified that this is equivalent 



\\CSG 

to (fT2"|) . Thus, (C, D) achieves the lower bound if and only 

if D = KC^ 1 and (fT2|) holds, since these conditions imply 

that (C,D)eQc,n( K )- 

Under the additional assumptions that K £ IZCi 

satisfies (JSJ, it will now be shown that there exists such 

(C, D) £ %2 x £2 with C outer. Since K satisfies (J6j) with 

M, N, Q,U,V £ WHoc it holds that 

log \K\ = log \MQ -U\- log \NQ + V\ 

By Theorem 17.17 in [27], \og\MQ - U\ £ d and 
log \NQ + V\ £ £1 and thus log \K \ £ C x . It follows from 
(fT3| and the boundedness of H, G yv and G zu on T that 



log 






-KH 



duj > — c 



and 






-KH 



£Ci 



Then by the factorization theorem in |30| there exists an 
outer function C £ H2 such that (fT2"j) holds. Also, D = 
KG- 1 £ £ 2 since 



\KC- 



\Kb HG zu G yv \^ 1 



KGy V 



HG ZU 



< 00. 



Remark 1: The spectral factorization gives some free- 
dom in the choice of (G, D) that attains the bound. For 
example, D instead of G could be chosen to be H2 and 
outer. That would result in having C £ £2- Considering 
more solutions than the one selected would require a 
slightly more complicated stability characterization, so 
this is not done. 

Remark 2: Optimal D will satisfy 



i n i2_ II"'-' HG zu G yv \\ 1 J G yv G* v 



c* c 



on T. 



It is interesting that the magnitudes of both G and D are 
directly proportional, on the unit circle, to the square root 
of the magnitude of K. In other words, the dynamics of a 
nominal controller K is "evenly" distributed on both sides 
of the communication channel. The static gain of G (and 
D) is tuned so that the SNR constraint is active. In the 
case when G 



y 



G zu , finding an optimal factorization 
amounts to performing a spectral factorization of \KH\ 
and tuning the static gain. If also H = 1 then the 
magnitudes of the frequency responses of G and D will 
then be proportional. 

C. Equivalent Convex Problem 

It will now be shown that a solution to the main problem 
can be obtained, with arbitrary accuracy, by solving a 
convex minimization problem in the Youla parameter. 

As discussed in the problem formulation, (G, D) should 
satisfy the SNR constraint (QJ and stabilize the system. 
The latter corresponds to T £ H2 or ([5]). Also, it was 
assumed that DC £ 1Z. Thus, the set of feasible (G, D) is 
given by 

®c,d = {(G, D) : DC £ n , ©, T £ H 2 } ■ 

Let M, N, U, V be determined by a coprime factorization 
of G yu and introduce 



A — M G zu G yv 

B = M 2 N- l VG zu G yv 

E = MNH 

F = (MV - 1)H 

L = G zv — al N G zu G yv . 



(14) 
(15) 
(16) 
(17) 
(18) 



It will now be shown that minimization of J(C, D) 
over ©en can be performed by minimizing the convex 
functional 

ip{Q) = \\L + AQ + B\\ 2 + — 2 , 

eH - \\EQ + F\\ 2 



over the convex set 



6o = 



[q-.Q 



£TIH 00 ,\\EQ 



F\\l< 



*} 



The Q £ Qq obtained from minimizing <p(Q) will be 
used to construct (C,D) £ &c,d- However, this will not 
be possible for Q for which the corresponding K has poles 
on T. For such Q a small perturbation can then be applied 
first. This will result in an increased cost, but this increase 
can be made arbitrarily small. That this is possible is 
established by the following lemma. 

Lemma 4: Suppose Q £ Qq and e > 0. Then there 
exists Q £ Qq such that 



and 



K = ^^£TIC U 

NQ + V 



ip{Q) < <p(Q) + e. 



(19) 



The proof of Lemma|4]is based on a perturbation argument 
and can be found in the Appendix. 

The main theorem of the paper can now be stated. 



Theorem 1: Suppose a > 0, that G yu = NM 



is 



a coprime factorization over TZHoo, that U, V £ TZHoo 
satisfy the Bezout identity VM + UN = 1 and that H £ 
IZWoo, G* ZU G ZU £ TZCoo and G yv G* yv £ TZCx, have no zeros 
on T. Then 



inf J(C, D) 

(C,D)£0 c ,n 



inf <p(Q). 



(20) 



Furthermore, suppose Q £ Qq, e > and let Q £ Qq be 
as in Lemma |4j Then there exists (G, D) such that the 
following conditions hold: 

. If MQ - U is not identically zero: (G, D) £ U 2 X £ 2 , 
where C is outer and 



K 



MQ-U 

NQ + V 



(21) 



|CT 



D 



KHG VU 



I- KG* 



K HGzuGyv 



(1 - KG yu f 
KG- 1 



IG* ZU G ZU Qn T 



GyvG yv 



(22) 
(23) 



. If MQ - U = 0: C = D = 0. 

If (G, D) satisfy these conditions, then (G, D) £ Qc,d and 

J(C,D) <ip(Q)+e. 

Proof: Consider (G, D) £ Q c ,d and define K = DC. 
Then (C,D) £ Qc,d(K) for this choice of K. Moreover, 
because T £ H2 it follows from Lemma Q] that K can 
be written using the Youla parametrization ([5]). Since the 
SNR constraint Q is satisfied by (G, D) it follows that 
K £ Qk, where 0# is defined by 



&K = <K : @, 



KHG yu 



1 - A'G„ 



< cr 



The inequality in this definition is strict because it was 
shown earlier that equality cannot hold. It has thus been 
proved that 

(C, D)£Q CD ^ (C, D) £ Qcd(K) for some K £ Q K . 

(24) 
A lower bound will now be determined for J(C, D). This 
will be accomplished through a series of inequalities and 
equalities, where each step will be explained afterwards. 



inf J(C, D) 

(c,fl)ee c , D 

> inf inf 

K£0k (C,D)£& c ,d(K) 



*~* zv 



-"- (-* ZU^ljV 



= inf 



> inf 
A'eej, 



G z 



(-*zv + 



■** ^- T zu^yv 



1-KG V 



-"■ ^ zu*~*yv 



1 — KG yu 



+ inf 

2 (C,D)eBc, D (K) 



DHG Z 



1 - ^G zu 

DHGzn 



\ — K. G yu 



KHGzuGyy 



I -KG, 

2 



(1 - AG yu ) 2 



KHG VU 



1 — KGy U 



= inf y>(Q) 

The first step follows from (f2"4")l and rewriting J(C, .D) 
in terms of K . In the second step, the first term has been 
moved out since it is constant in the inner minimization. 
The third step follows from application of Lemma [3] with 
2 



KG yU 



I- KG* 



>0, S 



1 



1-KG V 



£KH C 



2 ^ iiv^J/U 

Let A,B,E,F,L be given by (fl"4" ]) -([T5| ) . Application of 
the Youla parametrization and the Bezout identity then 
gives 



G Z v~\- 



■**■ ^ zu^yv 



G zv + 



1 



1 I G ZU Gy V G 



1 — KGy U \ 1 _ KGyu 

= G zv + (MNQ + MV - 1) GzuGyuMN- 1 

= AQ + B + G zv - MN- l G zu G yv 



and 



I\ 12 (sr zu (jr y V *~X ZU^yv^yu *^ *1 Kj y U 

(1 - KG yu Y = (1 - KGyuY 

= G zu G yv M 2 (Q + N~ X V) {MNQ + MV - 1) H 

= {AQ + B){EQ + F). 

The fourth step now follows from the definition of tp(Q). 
Now a suboptimal solution will be constructed. Suppose 
that Q £ Qq and e > and let Q £ Qq be as given by 
Lemma Hand define K £ KCi by d5J). Then K £ Q K 
and 



<p{Q) 



zv + 



1 - KG V 



KHG Z uG yV 



(1 - AG yu ) 2 



KHG y 



1 — AG™ 



If MQ — 17 = then A = 0, 



J(0,0) = ||G 2t ,||; = < /3 (Q)<^(Q) + 



and we are done. 

If, on the other hand, MQ — U is not identically zero 
then K is not identically zero. By Lemma [3J there then 
exists an outer C £ H2 and D £ £2 such that (|22l) and 
(l2"31) are satisfied. The lemma also says that such (C, D) 
satisfy 

2 



DiJG 2 



I -KG 



u » 



KHG zu G yv 



(1 



KGy U ) 



KHG,, 



and 



CG 



i)i' 



1 - KGyu 



<a l 



1 — KGyu 

KHG VU 



1 — KGyu 



D, C and if satisfy the conditions of Lemma [21 so T £ 
H2, which implies that (C,D) £ ©c,d- Moreover, 



J(C,D) 



G z 



KG ZU G. 



yv 



l-KGu 



DHG ZU 



1 — KGyu 



= v{Q) = <p(Q) 



Since e can be made arbitrarily small this shows that (j20l) 
holds and hence the proof is complete. ■ 

Remark 3: Theorem [T] shows that an e-suboptimal so- 
lution to the main problem can be found by minimizing 
tp(Q) over Qq. The obtained Q may have to be perturbed 
so that the resulting K has no poles on the unit circle. 
Then C is given by a spectral factorization and D is then 
obtained from C. 

A by-product of Theorem [1] is a necessary and sufficient 
criterion for the existence of a stabilizing LTI controller 
that satisfies the SNR constraint. 

Corollary 1: There exists (C,D) that stabilize the 
closed loop system of Fig. [5] subject to the SNR constraint 
((4]) if and only if there exists Q £ TZHoo such that 



\\{MNQ + MV-1)H\\\ <cr 2 

For the AWGN channel, we have that H 
condition can be written 



\MNQ - 



MV\\l < a 2 



1 



(25) 
1 and the 



(26) 



since MNQ + MV — 1 is strictly proper and thus orthog- 
onal to 1. 

Remark 4: Corollary [T] implies that the minimum SNR 
compatible with stabilization of a stochastically disturbed 
plant by an output feedback LTI controller with two 
degrees of freedom over an ACGN channel can be found 
by minimizing the left hand side of (|25|) over Q £ TZT-Loq. 

For the AWGN case, the analytical condition (J2J, 
presented in [6], is actually derived from a minimization 
of the left hand side of (f2l?|) . This means that the same 
condition is necessary and sufficient in the present problem 
setting as well, when the channel noise is white. This fact 
has been noted before in [11] . To elaborate, there is no 
plant disturbance in the setup of [6]. In that case, the SNR 
required for stabilizability will be the same regardless if the 
controller has one or two degrees of freedom. However, [31] 



considered the case when there is a plant disturbance and 
showed that the SNR required for stabilizability may then 
be larger than prescribed by ([2]). However, the controller 
in [3T] was assumed to only have one DOF (the encoder 
part was fixed to be a unity gain). Theorem 17 in [11J and 
this corollary shows that if the controller has two DOF, 
then ([2]) is again a necessary and sufficient criterion for 
stabilizability. 

For the ACGN case, however, this result is not identical 
to those in [12] and [14] since they assumed no plant 
disturbance and (effectively) controllers with one DOF. 

It will now be shown that the minimization of <p(Q) 
over Qq is actually a convex problem. To this end, define 
the functional 



p(a,e) = 



2tt 



a(ui) dui + 



*£ fl„ a(uj)e(uj)duj 

^ 2 + i-^OM^ 



with domain O p consisting of functions a and e that are 
continuous on [— 7T, 7r] and satisfy 



1 

2^ 



e(w) 2 dus < a 2 + 1. 



Lemma 5: The functional p(a, e) is convex. 
Proof: Take n > 2. The function 

f(x, y, v) = (x + yvf(x + yv) - v 2 , 

— x x + 2vx y + v (y y — 1) 

with domain Ux, y, v) : x, y £ R™, v £ R, y T y < l}, is 
convex in (x, y) for any v £ R. Thus, 



g(x,y) = max f fay, v) =x T x + 



{x T y) 
1-2/2/ 



with domain {(x,y) : x,y £ R™, y T y < l}, is convex in 
(cc, y) since it is the pointwise maximum of a set of convex 
functions [32]. Now, suppose (a, e) £ 9 p . Let 

wi = 0, uj k+ i - uj k = 2ir/n, fc = 1, ...,n — l 

a = [o(wi) o(w2) ••■ a(w„)] 

e = [e(wi) e(w 2 ) ■ ■ • e(w n )] . 

By definition of the integral, it holds that 

1 1 



lim 



~T - 

e e 



(a 2 + l)n (a 2 + 1) 2tt J_ n 



e(ui) du> < 1. 



So for large n, (a, (a 2 + 1) 1 / 2 e) / s/n belongs to the 
domain of g and 

p(a,e) = lim g -— , 



n y/{a 2 + l)n) 

Since the right hand side is convex in (d, e), and thus in 
(a, e), it follows that p(a, e) is convex. ■ 

Remark 5: Convexity of p(a, e) has been shown previ- 
ously in [33]. This proof is, however, substantially shorter. 

The convex functional p will be used in a relaxation of 
the minimization of <p(Q). A slight modification has to 



be done to tp(Q) in order to be able to compare the two 
functionals. For this purpose, define 

A(Q) = ||L||2+2Re(L,Ag + B>. 

Since A(Q) is affine in Q it doesn't affect the convexity 
of ip(Q). Define the functional 



<f>o{Q) = <P(Q) ~ A(Q) 

= \\AQ + B\\l + 



\(AQ + B){EQ + F)\\l 
a 2 - \\EQ + F\\l 



Lemma 6: Suppose Q £ Oq. Then (fo(Q) < 7 if and 
only if there exists (a, e) € O p such that p(a, e) < 7 and 



<H W ) > J G* ZU G zu G yv G* v 



M 2 Q 



M 2 V 



N 



e(u)>\EQ + F\, w€[-7r,7r]. 

Proof: Suppose first that <po(Q) < 7. Let 

M 2 V 



U G [— 7T,7r] 

(27) 
(28) 



a ( w ) — y G* zu G zu G yv G* yv 

e(uj) = \EQ + F\ 



M Z Q 



N 



and it follows that (a, e) S O p and p(a,e) = (fo(B,K). 
Conversely, suppose that (a, e) € P satisfy (j2"Tj) and (|2"5|) 
and that p(a, e) < 7. Then it follows from inspection of 
<Po{Q) and p(a,e) that tpo(Q) < p{a,e) < 7. ■ 

Convexity can now be proved. 

Theorem 2: The problem of minimizing y(Q) over 0q 
is convex. 

Proof: Suppose Qi, Q2 € @g- Then by Lemma[6]there 
exists (ai,ei) € O p and (02,62) € S P such that p(a\,ex) < 
yo(Qi) and /o(a2, 62) < yo(Q2)- It thus holds for < 9 < 1 
that 

9MQi) + (1 - ^)^o(Q 2 ) > 0p(a x , ei) + (1 - 0)p(a 2 , e 2 ) 
> p (6cn + (1 - 0)a 2 , 0ei + (1 - 0)e a ) 



The second inequality follows from Lemma [SJ The third 
inequality follows from Lemma [5] and that the constraints 
([27]) and ([28)) are convex. It is thus proved that <po(Q) is 
convex in Q. Then ip(Q) is convex since A(Q) is convex. 
It is finally noted that Qq is a convex set. ■ 

V. Numerical Solution 

By Lemma [6j minimizing tp(Q) over Qq is equivalent 
to minimizing /j(a, e) + A(Q) over p x Qq subject to 
(|27p and ([2"8")1 . This problem is infinite-dimensional, so the 
integrals are discretized for numerical solution. It will now 
be shown how the discretized problem can be posed as a 
semidefinite program. 

Let n > 2 and define {wfe} fc=0 , a and e as in the proof 
of Lemma [SJ Approximations of p(a, e) and A(Q) with n 



grid points are then given by 

1 



1:T: 



Pn(a,e) 



-a a + 



A„(Q) = ||L|| 



1 - ie T e 



V^ 



Re> tr (L*(AQ + B))\ z=e „ 



fe=i 



The accuracy of this approximation clearly depends on 
the number of grid points n. When implementing the 
minimization program, Q is parametrized using a finite 
basis representation. The accuracy of the approximation 
obviously depends on this representation as well. 

The denominator of p n (a,e) is positive for sufficiently 
large n and p n (a,e) + A„(Q) can be written as a Schur 
complement. It follows that p n (a,e) + A„(Q) < 7 if and 
only if 



ie T e -a 2 -I 
n 1~T~ 

—e a 



1 «r« 
-a e 



ia T a + A„(Q)-7 



^0, 



or, equivalcntly, 

'n(a 2 + 1) 








717- nA n (Q) 



[e a] T I [e a] y 0. 



Noting that the left hand side of the last inequality is also 
a Schur complement, it follows that it is equivalent to 



I e a 

e T n(a 2 + 1) 

a T wy — nA n (Q) 



>-0. 



(29) 



Let g k = ^G zu {e^)*G zu (e*"» ) G yv (e*"* ) G^ (e iw * ) * 
The constraints can then be approximated by 



icjfc ^2 I 



8(wi)>ftWe**) 



Q(e 2 ' 



F(e ia "=) 



N{e iu k) 



e(w fc ) > |£(e tWfc )Q(e Mfc ) + F(e^ fc )| 

1 - 

cr 2 + l > - Ve(w fe ) 2 . 



, k = 1 . . .n 

(30) 
l...n (31) 

(32) 



fc=i 



Minimizing 7 subject to ([29)) - ([32]) is a semidefinite pro- 
gram. 

A procedure for numerical solution will now be outlined. 

1) Determine N, M, U, V € IZHoo by a coprime factor- 
ization of G yu and calculate A, B, E, F and L. 

2) Choose n large, determine the grid points ujf., k = 
1 . . . n and solve the optimization problem of mini- 
mizing 7 subject to (|29l) - (|3"2")l . The transfer function 
Q is parametrized with a finite basis representation, 
for example as an FIR filter. If the problem is 
infeasible it could mean that a larger a 2 is needed to 
stabilize the plant. This can be checked analytically 
using the condition in [B]. If a 2 is sufficiently large 
according to this condition, the problem could still 
become infeasible if n is too small or Q is too coarsely 
parametrized. 

3) If NQ + V has zeros on the unit circle, determine a 
small perturbation Q of Q as outlined by Lemma [4] 

4) Determine K from (TT9")) . 
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5) Use a finite basis approximation A(uj) of CC* , for 
example the parametrization 



Control Performance 



N c 



A(u) = A + Y^ M 



^ — ki 



(33) 



fc=i 



and fit A{uS) to the right hand side of (|22p. for 
example by minimizing the mean squared deviation. 
Perform a spectral factorization of A(uS), choosing 
C as the stable and minimum phase spectral factor. 
7) Let D = KC~\ 



6) 



A. Example 

Consider the system in Fig. [3] A SISO plant is controlled 
over an AWGN channel. The SISO plant represents a 
special case where 



G{z) = 



G zv (z) 

Gy V (z) 



zu \Z ) 
yu\ z ) 



P(z) 
P(z) 



P(z) 
P(z) 



Let the plant be P(z) — l/(z(z — 2)). It has one unstable 
pole and a one-sample time delay. Using the stabilizability 
condition ([2]) . it is determined that stabilization is possible 
for a 2 > 12. (We have r\ = 0, since there are no non- 
minimum phase zeros, and r5 = 9, because of the location 
of the unstable pole and the relative degree, which is 2. 
For details, see [6]). 

A controller was determined for various values of a 2 , 
using the algorithm outlined above. The optimization was 
performed in Matlab, using the toolboxes Yalmip [34J and 
SeDuMi [35] . In the optimization program, n — 629 grid 
points were used and Q was parametrized as an FIR filter 
with length 20. The plant output variance is plotted in 
Fig. M for a number of different a 2 . It can be seen that the 
variance grows unbounded as a 2 approaches 12 and the 
feedback system comes closer to instability. This seems 
to be in agreement with the performance bound given in 

m 

VI. Conclusion 

This paper has considered a special class of decentral- 
ized control problems where the controller is split in two 
parts that are separated by a noisy communication channel 
with an SNR constraint. It has been shown that an optimal 
linear design can be obtained with arbitrary accuracy by 
solving a convex optimization problem and performing a 
spectral factorization. 
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Fig. 5. Control of a SISO plant over an AWGN channel. 
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Fig. 6. Minimum variance of the plant output y as a function of 
the SNR (or maximum allowed transmission power) a 2 , for the plant 
G = l/(z(z — 2)). The variance grows unbounded as a 2 approaches 
the lower limit for stabilization. 



The results in this paper can be viewed as a gener- 
alization of some results pertaining to a communication 
problem that can be obtained by considering the open- 
loop version of the control problem. 

As mentioned in Section HI the problem in this paper 
has previously been considered in the case with an AWGN 
channel with feedback |9j, where a similar result is ob- 
tained using a slightly different technique. A disadvantage 
with that result is, though, that it requires the controller 
to be over-parametrized with four degrees of freedom. The 
technique used in this paper has been applied to that case 
as well in [25], giving a solution that does not require an 
over-parametrization. 

Objects for further research include an extension to 
handle MIMO channels or plants with more than one 
controller input or measurement signal. Of course it would 
also be of interest to know if non-LTI controllers could 
provide better performance or require lower SNR levels 
for stabilization when the channel noise is colored. 

Appendix 

Lemma 7: Suppose that Y G M + is square and outer, 
X e Af+, and that Y^X e C p . Then Y^X G H p . 

Proof: Y -1 £ M + by Theorem 10 in [28]. It is easy to 
verify that the product of two J\f + functions is J\f + . The 
result follows from C p nM + = U. p [IS]- ■ 

Proof of Lemma|3J The proof is based on construction 
of Q through a perturbation of Q. Take Q e ©q and let 

MQ~U 



K 



NQ + V 



If K G TZC\ then let Q = Q and the construction is 
complete. Suppose instead that K has at least one pole 
on T. Since MQ — Ue IZHoo , z is a pole of K if and only 
if 

N{z)Q{z) + V{z) = 0. (34) 

Moreover, suppose that AM]) holds and that N(z) = 0. 
Then it follows from the Bezout identity that V(z) 7^ 0, 



which is a contradiction. Thus if NQ + V has a zero at z 
then N(z) ^ 0. 

Suppose now that NQ + V has a zero at zo € T and 
that zq ^ R (the case when zo € R is discussed later). Let 



Q = Q + A + Xiz \ A ,Ai e 



Then 



sg + F 



< a 2 + 1 if |Ao| + I Ail < 5 X for small 



enough S\. 

The coefficients Ao, Ai will be chosen so that the zero at 
zo is perturbed away from T. It must also be made sure 
that none of the other zeros can reach T under the same 
perturbation. For this reason, define the set of zeros not 
on the unit circle, 

n = {z : z <£ T,N(z)Q(z) + V(z) = 0}, 

and the smallest distance from that set to the unit circle, 



inf \zi 

zien,z 2 eT 



Z2 



where r > since Q has a finite number of elements. The 
location of the zeros of NQ + V depend continuously on 
(Ao, Ai). Thus, there exists S r > such that if |Ao| + |Ai| < 
S r then all zeros are displaced strictly less than r. 
Introduce the function 

X(z, Ao, Ai) = NQ + V = NQ + V + N(\ Q + Aiz" 1 ). 

Then 

dX' 



dct 



Re 
Im 
is non-zero at z 



dX 

9Ao" 
dX 

d\~o~ 



Re 



Irn 



dX 
9AT 



dct 



ReN ReA^ 1 
ImTV ImA^z" 1 



zq since N(zq) ^ and zq G T \ R. 
Then, by the implicit function theorem, there is a differ- 
entiable mapping z n- (Ao, Ai) defined in a neighborhood 
of zq, such that 

N{z)Q{z) + V{z) 
= N{z)Q{z) + V{z) + N(z)(X (z) + X^z' 1 ) = 0. 

This means that a new location z can be determined for 
the zero, and the mapping gives the corresponding Ao, Ai. 
Take e > 0. Since tp(Q) is continuous there exists Sq > 
such that 



Q-Q 



<6, 



Q 



<p(Q) - <p(Q) 



< £. 



Continuity of the mapping from z to (Ao, Ai) implies that 
there exists 5 Z > such that 

\z - z \ < S z =>■ |Ao(z)| + |Ai(z)| < mm{SQ,8x,5 r }. 

Now pick z ^ T such that | z — zq \ < S z and the mapping 
to Aq,Ai is defined. Then 



Q-Q 



< |A (z)| + |Ai(z)| < mm{S Q ,d\,S r }, 



which implies that 



EQ + F <<r 2 



\ip{Q) - <p{Q)\ < e, 



and that there are no new zeros on T. Since z ^ T it follows 
that NQ + V has at least one zero less than NQ + V on 
T. 

If zq is real, then define instead 

Q = Q + Xq, A gR 

and determine Ao analogously. Note, however, that the 
zero must be kept on the real axis. 

If Q is such that NQ + V has zeros on T, the procedure 
may be repeated again, with s appropriately chosen, until 
there are no such zeros. Thus, for every Q G Oq and e > 
it is possible to construct Q such that NQ + V has no 



zeros on T, \(p(Q) — <fi(Q)\ < £ and 



EQ + F 



<a 2 + l. 
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